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ABSTRACT 

We selected from the VLA FIRST survey a sample of 94 objects with starlike coun- 
terparts in the Sloan Digital Sky Survey, and with APM POSS-I colour O - E > 2, 
i.e. consistent with their being high-redshift quasars. 78 of the 94 candidates can be 
classified spectroscopically on the basis of either published data (mainly SDSS) or the 
observations presented here. The fractions of QSOs (51 out of 78) and redshift z > 3 
QSOs (23 out of 78, 29 per cent) are comparable to those found in other photometric 
searches for high-redshift QSOs. We confirm that selecting colour O — E > 2 ensures 
inclusion of all QSOs with 3.7 < z < 4.4. The fraction of 2 < z < 4.4 QSOs with broad 
absorption lines (BALs) is 27 ± 10 per cent (7/26) and the estimated BAL fraction 
for radio loud QSOs is at least as high as for optically selected QSOs (~ 13 per cent). 
Both the high BAL fraction and the high fraction of LoBALs among BALs (four to 
five out of 7) in our sample, compared to previous work, are likely due to the red colour 
selection O — E > 2. The space density of radio loud QSOs in the range 3.7 < z < 4.4 
(z = 4.0) with M AB (1450) < -26.6 and P > 10 25 7 W Hz" 1 is 1.7±0.6 Gpc~ 3 . Adopt- 
ing a radio-loud fraction of 13.4 ±3 per cent, this corresponds to p = 12.5 ±5.6 Gpc -3 , 
in substantial agreement with the cumulative luminosity function of SDSS QSOs in 
Fan et al. (2001a). We note the unusual flat-spectrum radio-luminous QSO FIRST 
1413+4505 (z = 3.11), which shows strong associated Lya absorption (rest-frame 
equivalent width ~ 40 A) and an extreme observed luminosity, L~2x 10 15 L Q . 

Key words: surveys - quasars: general - galaxies: high redshift - early Universe - 
radio continuum: galaxies 



1 INTRODUCTION 

This is the third of a series of papers presenting new samples 
of high-redshift radio QSOs selected by matching the FIRST 
catalogue of radio sources (Faint Images of the Radio Sky 
at Twenty-cm; Becker, White & Helfand 1995, White et al. 
1997) with red starlike objects from the APM (Automated 
Plate Measuring Facility) catalogue of the POSS-I survey 
(McMahon & Irwin 1992). 

Papers I and II (Benn et al. 2002, Holt et al. 2004) re- 
ported a sample of 18 z > 3.6 QSOs including the largest 
sample of z > 4 radio-selected QSOs then known. The search 
was carried out within a ~ 7030 deg 2 region using the con- 
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straints: i) E < 18.8 and starlike, ii) Sia GHz > 1 mjy, and 
iii) colour selection O — E > 3. This colour range includes an 
estimated 95 ± 1.5 per cent of E < 18.8 QSOs with redshift 
3.8 < z < 4.5 (Vigotti et al. 2003). On the basis of the 13 
QSOs with 2 ~ 3.8 — 4.5 and comparing to an equivalent 
sample at z ~ 2 drawn from the FIRST Bright Quasar Sur- 
vey of the north Galactic cap (FBQS-2, White et al. 2000), 
we showed (Vigotti et al. 2003) that the decline in space 
density of M A b (1450 A) < -26.9 QSOs (H a = 50 km s _1 
Mpc -1 , Q,m = 1 and Q.a = here and throughout the pa- 
per) was approximately a factor 2 between z ~ 2 and z ~ 4, 
significantly smaller than the value ~ 10 found for samples 
including lower luminosity objects (Fan et al. 2001a). 

The Sloan Digital Sky Survey Data Release 3 (SDSS, 
Stoughton et al. 2002; DR3 Abazajian et al. 2005) provides 
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Table 1. Right ascension and declination ranges of the four areas 
of sky included in the survey. 



RA (h) DEC (deg) 



7.933 


-> 8.667 


+30 -> 


+ 46 


8.733 


-> 9.600 


+38 -► 


+ 55 


9.667 


-> 15.667 


+42 -> 


+ 58 


15.733 


-> 16.667 


+30 -» 


+ 45 



moderately deep CCD imaging in five bands ugriz cover- 
ing ~ 5282 deg 2 . We present and discuss a new sample of 
high-redshift radio-selected QSO candidates in a 1378.5 deg 2 
area of overlap between FIRST and SDSS DR3 in the north 
Galactic cap. The selection criteria are: i) E < 19.1 and 
starlike in APM, ii) Si. 4 GHz > 1 mjy, iii) radio-optical sep- 
aration less than 1.5 arcsec, iv) colour O — E > 2 (includ- 
ing O non-detections) and v) starlike in SDSS. This new 
sample, with a wider colour range, has several advantages. 
Firstly, the SDSS photometric catalogue provides reliable 
morphological classification of the sources, allowing us to 
readily eliminate the galaxies classed as starlike in APM 
POSS-L Secondly, SDSS provides ~ 3 A-resolution spectra 
and spectroscopic classifications of many objects, particu- 
larly those selected as QSO candidates on the basis of their 
ugriz colours, or as counterparts of FIRST sources. We con- 
tinue to use the APM catalogue for colour selection since our 
previous work showed a high efficiency and completeness in 
the selection of z > 3.85 QSOs using O — E > 3 and with 
the new limit O — E > 2 we can check that no z > 3.7 QSOs 
have O - E < 3. 

The paper is structured as follows. In Section 2 we 
present the sample and the status of the spectroscopic clas- 
sification. Section 3 reports optical spectroscopy of part of 
the sample. The spectroscopic classification of the sample 
as QSOs, narrow emission line galaxies or stars is presented 
in Section 4.1. The distribution of optical magnitudes and 
O — E colours is discussed in Section 4.2. In Section 4.3 the 
sample is compared with previous radio-selected QSO sam- 
ples from the literature, in terms of the selection criteria 
and the resulting QSO redshift distribution. In Section 4.4 
we comment briefly on the spectra of seven QSOs exhibit- 
ing strong blueshifted broad absorption lines (BALs) and we 
analyse the fraction of BAL QSOs in the sample. Section 4.5 
is devoted to the peculiar QSO FIRST 1413+4505. In Sec- 
tion 4.6 we compute the absolute magnitudes, fc-corrections 
and radio luminosities of the 10 QSOs with z > 3.7, and we 
discuss the completeness of a sub-sample of seven of them. 
In Section 5 we use this sample to calculate the space density 
of QSOs. Section 6 summarizes our conclusions. 



2 SELECTION OF THE SAMPLE 

The sample was selected from the 1378.5-deg area defined 
in Table 1. This area includes most of the region covered 
by SDSS DR3 in the north Galactic cap, which is also cov- 
ered by the FIRST survey and by the APM catalogue of 
POSS-I. The FIRST survey includes 122463 sources in this 
area with S1.4 ghz (peak) > 1 mjy, of which 113 have APM 
E < 19.1, O — E > 2.0 and lie within 1.5 arcsec of starlike 
objects in SDSS DR3. Eighteen of these were undetected in 



APM POSS-I O but were detected in APS POSS-I O and 
had APS(0 — E) < 2, and were therefore removed from 
the sample. For the source FIRST 1340+5619 we found a 
large difference between the SDSS and APM magnitudes 
(r = 23.51 versus E = 18.99) and the source was eliminated 
from the sample after confirming with SDSS that the APM 
counterpart is a blend. The final sample thus includes 94 
candidate high-redshift radio QSOs. 

Of these 94 candidates, spectra were first obtained for 
seven in papers I and II, for six in the literature (found us- 
ing the NASA Extragalactic Database - NED) and for 41 
by SDSS DR3 (which also reobserved nine of the thirteen 
previously discovered). In Section 3 we present TNG optical 
spectroscopy of 13 (randomly selected) of the 40 remain- 
ing candidates, and we classify an SDSS DR3 source given 
spectral class 'unknown' in the Sloan survey. Subsequent to 
these observations, SDSS DR4 (2005 June 30) reported spec- 
troscopy of 10 of the remaining 27 candidates (and also of 
four of those observed here). In total, 78 of the 94 candi- 
dates (83 per cent) are now spectroscopically classified (see 
Section 4.1). 



3 TNG OPTICAL SPECTROSCOPY 

Spectra of 11 candidates (indicated in column 9 of Table 2) 
were obtained with the Telescopio Nazionale Galileo (TNG) 
on 2005 March 10 and 11 using the DOLORES (Device Opti- 
mized for LOw RESolution) spectrograph in long-slit mode. 
The LR-B grism was used, yielding a wavelength range 3000 
- 8800 A and dispersion 2.9 A pixel -1 . The detector was a 
thinned, back-illuminated Loral CCD with 15pim pixels. Ex- 
posure times were typically 900 s. Three spectrophotometric 
standard stars were observed in order to calibrate the instru- 
mental spectral response. The seeing was « 1.8 arcsec, and 
the width of the slit was set to 2 arcsec, yielding a spectral 
resolution of 23 A, as measured from sky lines. 

Standard data reduction was carried out using the 
IRAFt package. Arc-lamp exposures were used for the wave- 
length calibration, and yielded solutions with rms residu- 
als < 1.5 A . Small spectrum-to-spectrum wavelength shifts 
were corrected using sky lines in the science spectra. In most 
of the spectra, the two or three brightest sky emission lines 
were not completely removed in the sky subtraction. This 
did not affect spectral classification or redshift measurement, 
but for presentation in Figs. 1, 2 we manually cleaned some 
of the strongest residuals. 

Spectra of two more candidates (FIRST 1343+4305 and 
FIRST 1405+5155) were obtained with similar instrumental 
set-up on the nights of 2005 July 12 and 13. Two spectra 
per object were taken, shifting the object along the slit to 
minimize the effect of detector artefacts and remove the OH 
sky lines. For each object, one 2D spectrum was subtracted 
from the other, to remove the background, then wavelength 
calibrated, aligned in the spatial direction and coadded. Ex- 



< IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which is operated by Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Figure 1. TNG spectra of the four QSOs and the emission-line galaxy identified for the first time here, and the SDSS spectrum of source 
FIRST 1039+4931, which was classified 'unknown' in SDSS. The night sky spectrum (scaled) is shown near the zero level in the TNG 
spectra. The emission features are labelled by ion (wavelengths from Vanden Berk et al. 2001). 



posure times were 2 x 1800 s for 1343+4305 and 2 x 1200 s 
for 1405+5155. 

The 13 objects include eight QSOs, one z — 0.2705 
emission line galaxy (ELG), and four late-type stars. The 
redshift of each QSO was estimated as the average of the 
values measured from individual emission-line centroids (ex- 
cluding Lya, which is often affected by Lya forest absorp- 
tion). Six of the QSOs have 3.1 < z < 3.9, the other two 
have 2 = 0.897 and z = 0.943. 

The TNG spectra of the eight QSOs and the ELG are 
shown in Figs. 1 and 2. Fig. 1 also shows the SDSS DR3 spec- 
trum of FIRST 1039+4931, classified 'unknown' in SDSS, 
but which we identify as a QSO at z = 1.09. Four of the 
eight QSOs were also observed in SDSS DR4, published af- 



ter our TNG observations took place. The TNG spectra of 
these four objects are shown in Fig. 2. 



4 PROPERTIES OF THE SAMPLE 
4.1 Spectroscopic classification 

Table 2 lists the basic optical and radio properties of the 
94 objects in the sample, together with the spectroscopic 
classification and redshift, mostly from SDSS DR3 and DR4. 
The spectra were classified as one of: QSO (broad emission 
lines), emission- line galaxy (narrow emission lines), late- type 
star or early-type star. The sample of 94 objects includes 51 
QSOs (0.27 < z < 4.31), 25 stars, two ELGs (z = 0.27, 0.31) 
and 16 objects with no optical spectrum, i.e. 78 out of 94 (83 
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1110+4305 z = 3.B6 
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Figure 2. Spectra of four QSOs identified at the TNG and also observed by SDSS DR4. 




14 15 16 17 18 19 20 

E 

Figure 3. Distribution of corrected E magnitude. QSO identifica- 
tions are shown in black, stars and ELGs in grey and unclassified 
sources in white. 



per cent) are classified. Of the 78 objects with spectra, 38 per 
cent (30/78) are QSOs with z > 2, 29 per cent (23/78) are 
QSOs with z > 3 and 13 per cent (10/78) are QSOs with z > 
3.7, confirming the efficiency of the adopted selection criteria 
for identifying moderate- to high-z QSOs. The efficiencies of 
previous radio QSO surveys are compared in Section 4.3. 



4.2 Distribution of optical magnitudes and O — E 
colours 

In Table 2 and hereafter, unless otherwise indicated, the 
quoted POSS-I APM magnitudes are recalibrated with re- 
spect to APS (McMahon et al. 2002), and are corrected for 
Galactic reddening. This results in some of the sample hav- 
ing final O — E < 2 and E > 19.1 (the recalibrated E mag- 
nitudes are typically 0.3 ±0.3 mags fainter than the original 
APM values). The histogram of the E magnitudes is shown 
in Fig. 3. 

Fig. 4 shows the distribution O — E versus E. QSOs 
have a broader O — E colour range than stars. The ratio 
QSOs/stars increases from 21/12 in the range 17.1 < E < 
18.6 to 29/5 for E > 18.6 (see also Fig. 3). Assuming that 
the unclassified sources in Fig. 4 are either QSOs or stars, 
the distribution in O — E versus E suggests that ~ 60 per 
cent are QSOs. 

Fig. 5 (see also Fig. 6b) shows that the QSO colours 
redden significantly with redshift: the four measured values 
of O — E and three out of six lower limits imply O — E > 
3 for z > 3.7. This is due to the drop in intensity in the O 
band as it becomes more dominated by the region blueward 
of the Lyman limit, and the increase in intensity in the E 
band as Lya enters the band. 

An important conclusion from Fig. 5 is that the colour 
cut O-E > 2 includes most of the QSOs with 3.7 < z < 4.4, 
and this result is in agreement with previous studies. In 
Vigotti et al. (2003) we estimated a fraction of z > 3.8 QSOs 
with APM(0 — E) > 3 of 95 ± 1.5 per cent, using a sample 
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Figure 4. Colour-magnitude diagram for all 94 candidates in the sample. The dashed line corresponds to an assumed POSS-I blue plate 
limit O = 21.7 (referred to the APS calibration). Arrows show objects not detected in O band. 



Table 3. Selection criteria and QSO fractions for FBQS-2 and 
for the current sample 





FBQS-2 


This work 


Area (deg 2 ) 


2682 


1378.5 


E 


< 17.8 


< 19.1« 


O - E 


< 2 


>2« 


Radio-opt separation (arcsec) 


< 1-2 


< 1.5 


Starlikc morphology 


E or O 


E and SDSS 


Number of candidates 


1238 


94 


Number of spec, class. 


1130 


78 


Number of QSOs 


636 


51 


Number of QSOs 2 < z < 3 


53 


7 


Number of QSOs 3 < z < 3.7 


12 


13 


Number of QSOs z > 3.7 





10 



1: O and E magnitudes as taken directly from APM 
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Figure 5. O — E versus redshift for the QSOs in the sample. 



with E(APM) < 18.8. In FBQS-2 (White et al. 2000, O - 
E < 2 and E < 17.8) 90 per cent of the QSO candidates were 
spectroscopically classified, and these included 636 QSOs, 
none of which had z > 3.5, although 12 had 3 < z < 3.5. 
We are therefore confident that the number of QSOs we find 
with 3.7 < z < 4.4 is limited only by the radio and E band 
flux density limits, with no QSOs missed due to the colour 
selection. 



4.3 Comparison with other radio samples 

The survey reported here can be considered complementary 
to FBQS-2. Both surveys select FIRST radio sources with 
starlike optical counterparts, but the colour criteria are O — 
E < 2 for FBQS-2 (blue-excess objects) and O - E > 2 for 
our work. In FBQS-2 objects were selected with O — E < 2 to 
eliminate galaxies, since catalogues based on photographic 
plates provide poor discrimination between stellar and non- 
stellar objects. The survey selection criteria, and the QSO 
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Table 4. Well-defined samples of z > 3 radio-selected QSOs 



Reference, radio/opt limits Area (deg 2 ) Cands. Spec. Class. z > 3.0 z > 3.7 
0) (2) (3) (4) (5) (6) 

Hook et al. 1998 1600 73 50 6/50 12% 0/50 0% 

Green Bank S5 qhz > 25 mjy 
FIRST S 1A GHz cx\ 4 > -0.5 
APM POSS-I £ < 19.5, O - E > 1.2 

Snellen et al. 2001 6400 27 27 9/27 33% 5/27 19% 

Green Bank S5 GHz > 30 mjy 
NVSS S1.4 GHz af.4 > -0.35 
APM £ < 19, O - E > 2 



Hook et al. 2002 7265 BR 228 202 18/202 9% 12/202 6% 

Green Bank S 5 G Hz > 50 mjy (3637 BRI) 

NVSS Si. 4 GHz a? 4 > -0.5 

APM UKST R < 21; B - R > 1.5 

or_R-/>lorB-/>2 



Holt et al. 2004 7030 194 121 15/121 12% 15/121 12% 

FIRST Si. 4 GHz > 1 mjy 
APM E < 18.8, O - E > 3 



This work 1380 94 78 23/78 29% 10/78 13% 

FIRST Si. 4 GHz > 1 mjy 
APM E < 19.1, O - E > 2 

Col. (1): reference, optical and radio catalogues, and optical and radio flux-density and spectral-index (S„ oc u a ) limits. 
(2): survey area. (3,4,5,6): number of candidates, number of candidates with spectroscopic classification and the fraction of 
the latter at redshifts z > 3 and z > 3.7. 
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Figure 6. (a) O - E versus E for the QSOs in FBQS-2 (crosses) 
and in this work (circles for z < 3.7, squares for 2 > 3.7). The 
dashed line corresponds to the assumed POSS-I blue plate limit 
O = 21.7.(b) O — E versus z for the same QSOs. 



fractions as a function of redshift, are compared in Table 3. 
Figs. 6a and 6b show O — E versus E and O — E versus z for 
the QSOs in the two surveys. Whereas the redshifts covered 
by FBQS-2 range up to 3.4, our sample includes redshifts up 
to 4.3 (Fig. 6b, Table 3), and this is a consequence of both 
the colour cut and the fainter magnitudes. 

The efficiency of QSO selection in this work (51/78 = 
65 per cent) is slightly larger than in FBQS-2 (636/1130 = 



56 per cent). Although our selection criteria filter out most 
of the QSOs below z — 3.5, and these are covered by FBQS- 
2, our survey is better at rejecting galaxies, since we take 
the morphological information from SDSS whereas FBQS- 
2 used the less reliable APM discrimination. Although our 
sample has fewer sources, the efficiency for moderate to high 
redshift QSOs is high: 20/78 = 26 per cent for 2 < z < 3.7 
compared to 65/1130 = 6 per cent for FBQS-2. At z > 3.7 
the efficiency of our selection is 10/78 = 13 per cent. The 
maximum redshift in FBQS-2 is 2 = 3.4. 

Table 4 compares the sizes and search efficiencies of var- 
ious radio-selected high-redshift QSO samples. The success- 
rate for high-2 QSOs and the number of QSOs in this work 
are comparable to or larger than those of previous searches 
for radio-selected QSOs at much higher radio flux densities. 



4.4 The fraction of BAL QSOs 

Twenty-six of the QSOs with spectra have redshifts 2.0 < 
z < 4.4, so that Si IV 1397 A and C IV 1549 A are included 
in the observed wavelength range. Seven of these QSOs show 
strong broad absorption lines (BALs) with velocity widths 
above 2000 km s" 1 and reaching velocities above 3000 km 
s _1 , and were classified as BALs. The balnicity index (BI) of 
each BAL was computed following the prescription of Wey- 
mann et al. (1991) and is presented in Table 5. 

The fraction of BALs in our sample is therefore 7/26 
= 27 ± 10 per cent. For comparison, from the study of the 
BAL QSOs in FBQS-2 (Becker et al. 2000) we derive an ob- 
served fraction of BALs for 1.5 < z < 3 and 15 < E < 17.8 
of 14/134 = 10 ± 3 per cent (BALs with BI=0 excluded). 



© 2003 RAS, MNRAS 000, 1-?? 



A FIRST- APM-SDSS survey for high-redshift radio QSOs 7 




3000 7000 
Wavelength (A) 



000 7000 
Wavelength (A) 



Figure 7. SDSS spectra of four of the seven 2 < z < 4.4 BAL QSOs in the current sample. The three remaining spectra are presented 
in Fig. 1 of this paper, in Benn et al. (2002) and in Benn et al. (2005). 



In view of the redder colours found for BALs, compared to 
non-BALs, and the adopted colour selection, O — E < 2, 
Becker et al. conclude that the BAL fractions derived from 
their work should be considered as lower limits. Our result, 
27 ± 10 per cent, confirms a larger fraction of BALs in a 
red-selected sample (O — E > 2), and following their argu- 
ment this value should be considered as an upper limit to 
the actual fraction. We note however, that the distribution 
of balnicity indices for the FBQS-2 BALs reaches lower val- 
ues than the one in Table 5, which is more similar to the 
BI distribution of the BAL sample in Weyman et al. (1991), 
nowadays regarded as highly conservative. In particular we 
have 6/7 cases with BI > 4500 km s _1 (86 per cent), Wey- 
man et al. 24/40 (60 per cent) and Becker et al. 6/14 (43 
per cent). Because of our selection of the stronger BALs, in 
the comparison with FBQS-2 the fraction 27 ± 10 per cent 
should be regarded as a lower limit. Differences in the opti- 
cal magnitudes and redshift ranges, along with small number 
statistics could explain the differences in the BAL fractions 
between FBQS-2, > 10 ± 3 per cent, and our work, 27 ± 10 
per cent. 

Reichard et al. (2003a) obtained for a QSO sample se- 
lected from SDSS Early Data Release at 1.7 < z < 4.2 (re- 
gardless of radio emission) an observed BAL fraction 14 ± 1 
per cent. The BI distribution of these BALs reaches low 
values, with only 14/185 = 8 per cent of them having BI > 
5000 km s" 1 (Reichard et al. 2003b). Reichard et al. (2003a) 
confirm that BAL QSOs are redder than non-BALs, and es- 
timate a true BAL fraction of 13 ± 1 per cent, taking into 



Table 5. Balnicity indices of the z > 2 BAL QSOs 



FIRST 


2QSO 


BI 


BAL 


Spectrum 




name 




km s _1 


type 








0837+3641 


3.416 


5600 


LoBAL 


SDSS, 


Fig. 


7 


1134+4318 


2.184 


9800 


LoBAL? 


SDSS, 


Fig. 


7 


1219+4849 


3.201 


4800 


LoBAL 


TNG, 


Fig. 


1 


1408+5553 


2.446 


4500 




SDSS, 


Fig. 


7 


1459+4931 


2.370 


12850 


LoBAL 


SDSS, 


Fig. 


7 


1516+4309 


2.59 


6000 


LoBAL 


Benn et 


al. 2002 


1624+3758 


3.377 


2990 




Benn et 


al. 2005 



account colour-dependent selection effects in BALs and non- 
BALs. Trump et al. (2006) found from a larger sample of 
SDSS DR3 QSOs at 1.7 < z < 4.38 a slightly lower BAL 
fraction, 10.4 ± 0.2 per cent, with BI > 5000 km s _1 for 
162/1756 = 9 per cent of the BALs. For the comparison 
with our work and FBQS-2, the two last values should be 
regarded as upper limits, giving the low balnicity indices 
reached and the steeply rising distribution of BI towards low 
values (Reichard et al. 2003b). Taken together, the results 
from this work, Becker et al., Reichard et al. and Trump 
et al. suggest that the BAL fraction in radio-selected QSO 
samples is at least as high as that of optically-selected QSOs. 
A brief description of the seven BALs is given below: 

• FIRST 0837+3641 shows C IV and Si IV troughs start- 
ing at the QSO redshift, with a velocity width ~ 8500 km 



© 2003 RAS, MNRAS 000, 1-?? 



8 R. Carballo et al 



s _1 and sharp onset. The absorption bluewards of the ex- 
pected Al ill A1860 emission line suggests that this QSO is 
a low-ionization BAL (LoBAL). 

• FIRST 1134+4318 shows C IV and Si IV troughs start- 
ing at the QSO redshift. The C IV absorption presents two 
clearly distinct components. The weak absorption bluewards 
of the expected Al ill emission line suggests that this QSO 
could be a LoBAL. 

• FIRST 1219+4849 shows C IV and Si IV troughs, the 
former starting sharply at the emission redshift and extend- 
ing almost to the Si IV emission line. There is absorption 
bluewards of the Al ill emission line, so this is probably a 
LoBAL, although some contamination by sky bands cannot 
be excluded. 

• FIRST 1408+5553 shows C IV and Si IV troughs of 
width ~ 5500 km s~\ starting at the emission redshift. 

• FIRST 1459+4931 shows C IV and Si IV troughs with 
a sharp onset at the emission redshift. The C IV absorption 
extends to the Si IV emission line and clearly shows two 
components. The strong Al ill absorption confirms this QSO 
as a LoBAL. We identify the emission line at 7000 A as Fe ill 
UV 48 (rest frame wavelengths 2062.2, 2068.9 and 2079.65 
A, see Laor et al. 1997). 

• FIRST 1516+4309, with deep C IV and Si IV troughs 
was firstly reported as a BAL in Benn et al. (2002). The 
absorption bluewards of the expected Al ill suggests that 
this QSO is a LoBAL. 

• FIRST 1624+3758 is an unusual BAL (Benn et al. 
2005), with a clear detached C IV trough extending from 
-21000 to -29000 km s -1 . 

The fraction of LoBALs in our sample, four to five out 
of seven or 55-70 per cent, is unusually high. For comparison, 
the observed fractions obtained from the FBQS-2 and from 
the SDSS BALs in Reichard et al. (2003a) are 28 ± 14 per 
cent (four out of fourteen) and 13 + 3 per cent (24/181, 1.7 < 
z < 3.9) respectively. The higher fraction in our sample is 
consistent with Becker et al. result that BALs in general 
and the LoBALs in particular are over-abundant among the 
reddest QSOs. Moreover, Reichard et al. (2003b) found that 
LoBALs have stronger CIV absorption troughs (i.e., stronger 
BI) than HiBALs, and the balnicity indices in our sample 
are in fact higher than the typical values in FBQS-2 and in 
Reichard et al.(2003a) SDSS sample. 

4.5 The peculiar QSO FIRST 1413+4505 

FIRST 1413+4505 (Fig. 2, z = 3.107) has an unusual spec- 
trum, with a broad strong absorption starting at the ex- 
pected position of the Lye* emission line and extending up 
to —18500 km s _1 bluewards. We identify the feature as H 
I on the basis of the location of the line and the detection of 
Ly/3 absorption at this redshift. The broad absorption is not 
detected in metal lines. The Lye* absorption feature has two 
components. The first one, at z a b s = 3.06, has a rest-frame 
EW ~ 38A and associated Ly/3 absorption with EW ~ 5A. 
The second component, at z abs — 2.93, has EW ~ 7A and 
the corresponding Ly/3 absorption is difficult to identify due 
to the high noise in this part of the spectrum. The EW of 
the first Lye* absorber is consistent with a single damped 
Lye* system with high H I column density, i.e. > 10 21 cm~ 2 
(Lanzetta et al. 1991), at the QSO position. However, the 



Lye*/Ly/3 ratio suggests an optically thin absorber, there- 
fore the system could be a velocity-resolved BAL-like wind. 
A higher resolution spectrum is needed to determine the 
nature of this absorption system. 

The spectral energy distribution of the source can be 
obtained from NED. The FIRST integrated flux density 
at 1.4 GHz is 140.3 ± 0.1 mjy, and the flux density at 
4.85 GHz is 125 ± 15 mjy (Becker, White & Edwards 
1991). The source is unresolved in FIRST, < 1 arcsec. 
The SDSS magnitudes, corrected for Galactic extinction, 
are u = 22.32 ± 0.23, g = 20.15 ± 0.04, r = 19.27 ± 0.02, 
i = 19.09 ± 0.02 and z = 19.09 ± 0.05. The source is not 
detected in 2MASS and the magnitude limits for 99% com- 
pleteness are: J = 15.8, H = 15.1 and K = 14.3 (see 2MASS 
Web sitet). FIRST 1413+4505 is highly luminous in the ra- 
dio, Pi. 4 GHz = 2.2 x 10 27 W Hz -1 , and has a flat radio 
spectral index, a = —0.09. From the E magnitude and us- 
ing the TNG spectrum for the fc-correction we obtained Mab 
(1450 A) = -25.7 (see equations in Vigotti et al. 2003). 

At the position of the QSO, mid-infrared emission is 
detected via ID addition of IRAS Scans (SCANPI) at 12^m 
and 25/im. The automatic procedure obtained via NED gives 
an 11a detection at 12/im in the addition modes 'detec- 
tor weighted' and 'mean', 15cr detection for 'rms weighted' 
and non-detection for the 'median' mode. The correspond- 
ing peak flux densities in Jy are 0.22 ± 0.02, 0.21 ± 0.02 and 
0.31 ± 0.02, respectively. At 25/im coaddition gives 3.6<r de- 
tection for 'detector weighted' and 'mean' modes, 5er for 'rms 
weighted' and non-detection with the 'median' mode. The 
corresponding peak flux densities are 0.08 + 0.02, 0.08 + 0.02 
and 0.11+0.02. 

Mid-infrared detection at z > 2 is rare and we checked 
if the detection with SCANPI could arise due to confusion 
with one of the two bright stars 52 arcsec NE and 48 arcsec 
SW (binary) of the QSO, but neither of them was detected 
at 12pm or 25pm. 

Assuming that the mid-infrared emission is physically 
associated with the QSO, and not a chance coincidence, the 
luminosity of the QSO in the rest frame wavelength range 
2200 A - 6.1 pm (from 2-band to 25 /jm in the observer 
frame) would be 1.8 x 10 15 L Q , i.e. FIRST 1413+4505 would 
be one of the most luminous objects known. To our knowl- 
edge there are five other systems with such extreme lumi- 
nosities, and although they were discovered in different ways, 
and have different SEDs and absorption/emission-line prop- 
erties, the extreme luminosity was in all cases ascribed to 
flux magnification due to gravitational lensing. These ob- 
jects are the Seyfert 2 like object IRAS FSC 10214+4724 
(Rowan-Robinson et al. 1991), the BAL QSO H1413+117 
(Magain et al. 1988), the sub-mm galaxy SMM 02399- 
0136 (Ivison et al. 1998), the optically bright QSO APM 
08279+5255 (Irwin et al. 1998) and the Lyman break galaxy 
MS 1512-cB58 (Yee et al. 1996). 

4.6 The z > 3.7 QSO sample 

Fig. 5 suggests that no QSOs with z > 3.7 will have been 
missed by the colour selection O — E > 2. The properties 
of the 10 z > 3.7 QSOs are given in Table 6. Adopting 

t http: / / www.ipac.caltech.edu/2mass 
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Table 6. Absolute magnitudes, radio luminosities and contributions to space density, of the z > 3.7 QSOs 



RA 


DEC 


2 


E 


51.4 


Mab 


log Pi. 4 




^max.o 


■^iiiax, r 


2 a 


p 


J2000 


J2000 






mjy 


1450 A 


W Hz" 1 










Gpc- 3 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


09 41 19.44 


51 19 33.0 


3.850 


18.02 


2.49 


-27.67 


26.00 


3.70 


4.40 


4.40 


3.70-4.40 


0.145 


10 57 56.28 


45 55 53.1 


4.137 


16.87 


1.10 


-28.37 


25.70 


3.70 


4.40 


4.38 


3.70-4.38 


0.149 


11 02 01.91 


53 39 12.7 


4.297 


19.64 


4.45 


-25.84 


26.33 












11 10 55.22 


43 05 10.1 


3.862 


18.08 


1.21 


-27.41 


25.69 


3.70 


4.40 


4.32 


3.70-4.32 


0.162 


12 40 54.92 


54 36 52.2 


3.938 


19.49 


15.09 


-25.85 


26.80 












13 09 40.70 


57 33 09.9 


4.274 


18.36 


11.33 


-26.59 


26.73 


3.85 


4.40 


4.40 


3.85-4.40 


0.187 


15 06 43.81 


53 31 34.4 


3.790 


17.89 


14.63 


-27.57 


26.76 


3.70 


4.40 


4.40 


3.70-4.40 


0.145 


15 10 02.93 


57 02 43.4 


4.309 


19.23 


254.97 


-26.53 


28.09 












16 19 33.65 


30 21 15.1 


3.794 


18.70 


3.88 


-26.96 


26.18 


3.70 


4.40 


4.40 


3.70-4.40 


0.145 


16 39 50.52 


43 40 03.7 


3.990 


17.68 


25.23 


-27.92 


27.03 


3.70 


4.40 


4.40 


3.70-4.40 


0.145 



Cols. (1, 2, 3, 4, 5): similar to cols. (1, 2, 11, 6, 4) in Table 2. (6): absolute magnitude Mab at rest-frame 1450 A. (7): log 
radio luminosity at rest-frame 1.4 GHz. (8, 9, 10, 11): rcdshift limits in the optical (2 m i njOI 2 m ax,o), radio (2max,r), or both 
(2 a ), within the range 3.7 < z < 4.4. (12): contribution to the space density. 



the magnitude limit E < 19.1, leaves 7 QSOs in the range 
3.7 < z < 4.4. 

The absolute magnitudes at rest-frame 1450-A, 
Mab (1450 A), of the z > 3.7 QSOs were calculated from 
the extinction-corrected E magnitudes using the same pro- 
cedure as in Vigotti et al. (2003), except for a better estima- 
tion of the fe-corrections, obtained here from the individual 
spectra rather than using an average. We used the optical 
spectra from this work (FIRST 1110+4305), Benn et al. 2002 
(FIRST 0941+5119) or from SDSS (remaining QSOs). 

Fig. 8a shows m-AB [1450(1 + z)] — E, or fc-correction, 
for the redshift range 3.7 - 4.4, using the spectra of the 10 
QSOs. Fig. 8b shows the mean and standard deviation of 
m AB [1450(l + z)]-E for the seven QSOs with E < 19.1. 

The radio luminosities at rest-frame frequency of 1.4 
GHz were calculated assuming a spectral index a = —0.3 
(£„ oc v a ), which is the median spectral index obtained for 
13 z > 3.6 FIRST- APM QSOs in Holt et al. (2004). The 
radio luminosities range from Pi. 4 ghz = 10 25 ' 7 to 10 28 ' 1 W 
Hz -1 , so all ten QSOs are 'radio loud' QSOs, adopting the 
criterion Pi. 4 GHz > 10 25 5 W Hz" 1 (Gregg et al. 1996). 

4.6.1 Completeness of the z > 3.7, E < 19.1 QSO 
sub-sample 

There are several sources of incompleteness, summarised be- 
low. 

(1) Only 78 of the 94 candidates were spectroscopically 
classified, giving a fraction of 83 per cent. 

(2) The APM completeness for E < 19.1 was estimated 
as the fraction of SDSS r < 19.8, z > 2 QSOs in the SDSS 
DR3 Quasar Catalog (Schneider et al. 2005) detected and 
starlike in APM E. The adopted r band limit was obtained 
from the average magnitude difference r — E = 0.69 (stan- 
dard deviation 0.34) of the 10 high-2 QSOs in the sample. 
Average r — E and standard deviation for the 94 candidates 
are 0.55 and 0.36 respectively. From the 2060 DR3 QSOs in 
the surveyed area of this work, 2002 were detected in APM 
and 1793 of these were starlike, giving a fraction 87 ± 2 per 
cent. 

(3) The completeness of the FIRST catalogue for the 
range 1 - 1.5 mjy is ~ 70 per cent (Prandoni et al. 2001). 




3.7 3.8 3.9 



4 4.1 4.2 
redshift 



1.5 



i 0.5 



-0.5 




7 QSOs E < 19.1 

mean and standard deviation 



3.7 



3.9 



4 4.1 
redshift 



4.2 



4.3 



4.4 



Figure 8. (a) fc-correction between the AB magnitude at 1450 
(I + 2) A and through the E band, versus rcdshift, obtained from 
the spectra of the 2 > 3.7 QSOs. Solid lines correspond to the 7 
QSOs with E < 19.1. (b) Average fc-correction for these 7 QSOs 
(solid line) and standard deviation (dashed lines). 



For the sample of seven QSOs, with radio flux densities: 
1.1, 1.2, 2.5, 3.88, 11.33, 14.6 and 25.2 mjy, we estimate a 
completeness of 92 ± 5 per cent, assuming Poisson errors. 

(4), (5) In Vigotti et al. (2003) we found for a very sim- 
ilar sample completeness 98 ± 1 per cent due to QSOs prob- 
ably missed because of extended radio emission, and 99 per 
cent completeness due to QSOs that may exceed the radio- 
optical separation limit of 1.5 arcsec. 

(6) From our work and from the literature (see Figs. 5, 
6 and Section 4.2) we estimate 100 per cent completeness 
for the colour cut O — E > 2 in the redshift range 3.7 - 4.4. 

(7) The completeness of the SDSS photometric survey 
quoted in the project web site, computed by comparing the 
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number of objects found by SDSS to the number found by 
the COMBO survey (Classifying Objects by Medium-Band 
Observations^), is 100 per cent for point sources with SDSS 
r < 21, which is above the maximum r = 20.6 in our sample. 

The combined completeness due to the above seven fac- 
tors is 64 ± 5 per cent. 



5 QSO SPACE DENSITY AT HIGH-REDSHIFT 

The observed space density of QSOs in the redshift range 
3.7 < z < 4.4, over the 1378.5 deg 2 survey area, was cal- 
culated using the sub-sample of seven QSOs with E < 19.1 
and the 1/V a estimator (see e.g. Avni & Bahcall 1980). The 
space density contributed by each QSO was computed as 
the inverse of the available volume, using the constraints 
E < 19.1, Sia GHz > 1 mJy and 3.7 <z < 4.4. The redshift 
limits and the space density contribution of each QSO are 
listed in Table 6. The sum of these contributions yields a 
space density (for M AB (1450) < -26.6) of 1.1 ± 0.4 Gpc~ 3 , 
assuming Poisson errors. The mean redshift of the QSOs is 
4.0, the mean absolute magnitude Mab(1450) = —27.5 and 
the mean log radio luminosity Pi. 4 ghz = 10 26 ' 3 W Hz -1 . 
Correcting for the 64 ± 5 per cent completeness, the space 
density is 1.7 ± 0.6 Gpc -3 . Assuming the radio-loud frac- 
tion 13.4 ± 3 per cent for P1.4 GHz > 10 25 ' 7 W Hz -1 from 
Vigotti et al. (2003), the space density of all QSOs with 
3.7 < 2 < 4.4 and M AB (1450) < -26.6 is 12.5 ± 5.6 Gpc~ 3 . 

5.1 Comparison with other surveys at 2 > 3.7 

In Vigotti et al. (2003) we computed p[M A b(1450) < 
-26.9] = 7.4 ± 2.6 Gpc~ 3 , at 3.85 < 2 < 4.45, assuming 
the same radio-loud fraction as above. The average redshift 
(4.2), optical luminosity (Mab(1450) = -27.7) and radio lu- 
minosity (Pi. 4 ghz = 10 26 ' 6 W Hz -1 ) were all slightly larger 
than for the current sample. Although the space density de- 
rived from the new sample has a larger statistical error, this 
work extends the computation of the space density to lower 
optical luminosities and benefits from better fc-corrections. 

Four of the seven z > 3.7 QSOs in the current work 
were also included in Vigotti et al. The Mab(1450) values 
obtained in the two analyses agree within 0.15 magnitudes, 
except for FIRST 1309+5733, with M AB (1450) = -26.6 
compared to —27.1 in Vigotti et al. The reason for the dis- 
crepancy is the high fc-correction for this source (the highest 
among E < 19.1 QSOs in Fig. 8). 

Fan et al. (2001a,b) obtained the space density of QSOs 
at high redshift using an SDSS-selected sample of 39 QSOs 
with 3.6 < z < 5.0 and -27.75 < M AS (1450) < -25.5. 
Seven of these QSOs have redshifts and absolute magnitudes 
within the range used here (2 = 3.7 - 4.4, Mab < —26.6) 
and in this sense our sample, with seven radio loud QSOs, 
has comparable quality in terms of the number of objects. 

Fig. 9 shows the 1/V a estimators of the cumulative lu- 
minosity function obtained in this work and in Vigotti et al., 
along with the maximum likelihood solution by Fan et al. 
For 2 = 4.0 Fan et al. give p[M AB (1450) < -26.6] = 8.0±4.6 

§ http: / /www. mpia-hd.mpg.de/COMBO 
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Figure 9. Cumulative luminosity function derived from the l/V a 
estimator in this work (circle, 3.7 < z < 4.4) and in Vigotti et 
al. (square, 3.8 < z < 4.5). Solid and dashed lines: maximum 
likelihood solution and standard deviation from Fan et al. (3.6 < 
z < 5.0). 

Gpc~ 3 , compared to 12.5 ± 5.6 Gpc~ 3 in this work. These 
two values are consistent, within the errors. 



6 CONCLUSIONS 

A sample of 94 radio-emitting QSO candidates has been de- 
rived by cross-correlating the FIRST survey, APM POSS-I 
and the SDSS photometric survey, selecting objects with 
E < 19.1, O - E > 2, and starlike in SDSS. 78 of the 94 
sources (83 per cent) are spectroscopically classified, mainly 
from SDSS, but also from the literature and from TNG spec- 
troscopy presented here (13 sources). The classified sources 
include 51 QSOs, with redshifts 0.27 < z < 4.31. Our main 
results are as follows: 

(i) The efficiency of selection of high-redshift QSOs is 
23/78 = 29 per cent for z > 3, and 10/78 = 13 per cent 
for 2 > 3.7, comparable to that for previous searches for 
high-redshift QSOs at much higher radio flux densities. 

(ii) We found no 3.7 < z < 4.4 QSOs with 2 < O - E < 
3, supporting our assumption that no z > 3.7 QSOs are 
missed by selecting O — E > 2. 

(iii) The BAL fraction for the QSOs with 2 < z < 4.4 
in our sample is ~ 27 ± 10 per cent (7/26), larger than 
found for FBQS-2 (Becker et al. 2000) and for the optically- 
selected SDSS samples in Reichard et al. (2003) and Trump 
et al. (2006), despite these three works being more complete 
for lower balnicity index. Our sample also has an unusually 
high fraction of LoBALs, four to five out of seven, compared 
to the above samples. The most likely explanation for the 
high frequency of BALs and LoBALs in our sample is the 
red colour selection, O — E > 2. A comparison of the three 
studies suggests that the BAL fraction in radio-selected QSO 
samples is at least as high as in optically-selected ones. 

(iv) We note the unusual QSO FIRST 1413+4505 (z = 
3.11), whose optical spectrum reveals a strong associated 
Lye* absorber (damped or BAL-like) that completely re- 
moves the Lya emission. The source is highly luminous in 
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the radio, has a flat radio spectral index, and is detected in 
IRAS 12 and 25 fim, yielding a bolometric luminosity in 
the range from SDSS z band to 25 fim L ~ 1.8 x 10 15 L©, 
which places it amongst the most luminous objects known. 

(v) Using a sub-sample of seven QSOs with 3.7 < 
z < 4.4 and E < 19.1, we derive a space density 
of radio loud QSOs, p[M AB (1450) < -26.6, P 1A gh z > 
10 25 ' 7 (W Hz" 1 )] = 1.7 ±0.6 Gpc~ 3 . Assuming a radio-loud 
fraction 13.4 ± 3 per cent, the space density for all QSOs is 
p[M A b(1450) < -26.6] = 12.5 ± 5.6 Gpc~ 3 . This result is 
in good agreement with the best-fitting model in Fan et al. 
(2001a, b), considering the statistical errors. 
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Table 2. Radio/optical properties and spectroscopic classification of the 94 candidates 
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Table 2. Continued 
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Table 2. Continued 



RA DEC FIRST-SDSS Si. 4 E(B - V) E O - E Spectroscopic classification Notes 

J2000 J2000 arcsec mjy SDSS Other Class z 
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Paper I 
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Cols. (1,2): SDSS coordinates. (3): radio-optical (FIRST-SDSS) separation. (4): FIRST 1.4 GHz flux density (integrated for 51.4 GHz > 10 mjy, otherwise peak, see fig. 4 
in Ivezic et al. 2002). (5): E(B - V) Galactic extinction (Schlegel, Finkbeiner & Davis 1998). (6,7): E magnitude and O - E colour taken from McMahon et al. (2002) 
FIRST-optical identification catalogue (APM magnitudes recalibrated with respect to APS, for more accurate photometry), and corrected for Galactic reddening assuming 
A(E) = 2.7E(B — V) and A(0) = 4AE(B — V) (see White et al. 2000); (8,9,10,11): source of spectrum, spectroscopic classification, and redshift. The references from the 
literature arc as follows: FIRST 1011+4239 and FIRST 1502+4837, Snellen et al. 2001; FIRST 1039+4705, Kopylov, Artyukh & Kopylova 2004; FIRST 1057+4555, Peroux 
ct al. 2001; FIRST 1510+5702, Storrie-Lombardi et al. 1996; and FIRST 1556+3517, Najita, Dey & Brotherton 2000. (12): Notes as follows: a = Not a QSO in the "DR3 
Quasar Catalog", probably because it doesn't meet the requirement Mj > —22 (H = 70 km s Mpc~ , !)m = 0.3, and Q\ = 0.7), b = Classified in this work from the 
presence of CIII] and Mgll emission lines in the SDSS DR3 spectrum, c = Spectrum not published, d = The source is not included in McMahon et al. (2002), therefore the 
APM calibration was adopted, e = spectrum not available on the SDSS website, f = See Section 4.5, g = Classified as a galaxy in 2MASS (The Two Micron All Sky Survey, 
Klcinmann et al. 1994) 



